Few-Body Systems 0, 1-11 (2008) P 

bystems 



© by Springer- Vcrlag 2008 
Printed in Austria 



Electromagnetic properties of diquarks 

o 

O ■ p. Maris* 

(N 

^ ' Dept. of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA 15260 

cn 
(N 

(N 



Abstract. Diquark correlations play an important role in hadron physics. 
^ ' The properties of diquarks can be obtained from the corresponding bound 

OO . state equation. Using a model for the effective quark-quark interaction that 

I has proved successful in the light meson sector, we solve the scalar diquark 

Q^ ' Bethe-Salpeter equations and use the obtained Bethe-Salpeter amplitudes to 

Cp . compute the diquarks' electromagnetic form factors. The scalar ud diquark 

(<~^ I charge radius is about 8% larger than the pion charge radius, indicating that 

these diquarks are somewhat larger in size than the corresponding mesons. 
We also provide analytic fits for the form factor over a moderate range in 
^ , Q^, which may be useful, for example, in building quark-diquark models of 

^ ' nucleons. 

c ■ 
> 

^ I 1 Diquarks 



Diquarks are coloured states of two quarks. Due to their colour charge, they 
are expected to be confined, and cannot be observed directly. Despite being 
confined, they can play a role in the dynamics of baryons: two quarks in a colour 
anti-triplet configuration can couple with a third quark to form a colour singlet 
baryon. Indeed quark-diquark models have been used quite successfully in various 
models for baryons |^ HI I3 . More recently, they also have been suggested as the 
constituents of pentaquarks: one could imagine a pentaquark as a bound state 
of two diquarks and an antiquark 4 , or as a bound state of a diquark and a 
triquark 5^. 

For calculations of the baryon spectrum within a quark-diquark model, the 
most important property of the diquarks is their mass. There are a number of 
different calculations of these masses available, both from lattice QCD jOj and 
using the set of Dyson-Schwinger equations [DSEs] [Tj |H1 El • These calculations 
seem to agree with each other to within about 20%. For calculations of elec- 
tromagnetic properties of nucleons (and other baryons), one also needs to know 
how photons couple to diquarks. Since they are not observable, we cannot obtain 
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this information from experiments (at least not directly), and we have to rely on 
theoretical calculations of the electromagnetic properties of diquark correlations. 

Here we discuss the electromagnetic form factors of the scalar diquarks, using 
the DSE approach |1U1 1111 1121 113j . The model we use for the interaction has 
been shown to give an efficient description of the light pseudoscalar and vector 
mesons |141 115j . In particular the pion and kaon electromagnetic form factors are 
in excellent agreement with the available experimental data [El • Our calculation 
can serve as a guide to further improvements of quark-diquark models of nucleons 
and other hadrons. 

2 Dyson Schwinger Approach 

The DSE for the renormalized quark propagator in Euclidean space is 

rA 

S{p)-' = iZ2^ + Z4mg{fi) + Zi g^D^,{k)>^^^S{q)rl{q,p) , (1) 

Jq 

where D^^{k = p — q) is the dressed-gluon propagator, ri{q,p) the dressed- 
quark-gluon vertex with colour-octet index i = 1,...,8, and Z2 and Z4 are 
the quark wave-function and mass renormalization constants. The notation 

= J d^g/(27r)^ stands for a translationally invariant regularization of the in- 
tegral, with A being the regularization mass-scale. The regularization is removed 
at the end of all calculations, by taking the limit yl ^ 00. We use the Euclidean 
metric where {7^,7;^} = 25^jy, jjt = 7^ and a ■ b = X]i^=i o-ih- The most general 
solution of Eq. has the form S{p)~^ = i'^A{p'^)+ B{p'^) and is renormalized 
at spacelike /i^ according to A{fj?) = 1 and -B(/i^) = rnq{^) with ?Tig(/x) being the 
running current-quark mass. 

For practical calculation, we have to make a truncation of the set of DSEs |1U1 
[n^lE]- Here we adopt the rainbow truncation of the quark DSE, given by 

Z,g''D,,{k)n{q,p) - , (2) 

where D^^^'^{k = p — q) is the free gluon propagator in Landau gauge and Qik"^) 
a phenomenological effective interaction |14| II 5j . 

2.1 Diquarks 

Diquarks are coloured quark-quark correlations. Two quarks can be coupled in 
either a colour sextet or a colour anti-triplet. Single gluon exchange leads to 
an (effective) interaction that is attractive for diquarks in a colour anti-triplet 
configuration, but the interaction is repulsive in the colour sextet channel [Jj. 
Furthermore, it is the diquark in a colour anti-triplet that can couple with a 
quark to form a colour-singlet baryon. Thus we will only consider the diquark 
correlations in a colour anti-triplet configuration. 

Using the antisymmetric tensor Cabc as a representation of the anti-triplet, the 
diquarks are described by solutions of the homogeneous Bethe-Salpeter equation 
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[BSE] 

{rD{p+,-p.)cY^\a,, = -e,,, j\'^2l\p^,q^-,-p.,-q.)C'"' C^^' 

S^^'{q+) [rDiq+, -q-)cY'' Siq.f" . (3) 

where = p + P/2 and —p~ = —p + P/2 are both outgoing quark momenta 
(and similarly for qj^ and —q~.), and C is the charge conjugation matrix C, which 
satisfies = —1 and Cj^C = {'^^Y' ■ 

In ladder truncation this equation reduces to 

rn{p+,-p-)C = j''g{k')D'^'^{k)^^s{q+)(rD{q+,-q-)c)s{q^)^,. 

(4) 

This equation resembles the bound state BSE for a meson, the only difference 
being that the effective interaction in the diquark BSE is reduced by a factor of 
two compared to that in the meson BSE in ladder truncation. 

The allowed flavour structure of the diquarks follows from the Pauli principle. 
The total wave function has to be antisymmetric in all its labels, and we know 
already that they are antisymmetric in their colour indices. For three flavours 
that implies that the scalar diquarks have to be in a flavour anti-triplet. On 
the other hand, axial-vector diquarks are in a flavour sextet. Convenient flavour 
matrices for the anti-triplet states are the Gell-Man matrices A2, — A5, and A7. 
Explicitly, that means that we consider only ud, us, and ds scalar diquarks here. 

Different types of diquarks are characterized by different Dirac structures. 
Since the intrinsic parity of a quark-quark pair is opposite to that of a quark- 
antiquark pair, the scalar diquark Bethe-Salpeter amplitude [BSA], or rather 
rD,scaiiQ+, —<1-)C , has the general decomposition 

rD,scAQ+,-q-)C = l^[iE{q\q-P)+ If>F{q\q-P) 

+ ^ G{q\q ■ P) + cj^, P^q, H{q\ q ■ P)] , (5) 

which is the same form as a pseudoscalar meson BSA. 



2.2 Electromagnetic interactions 



The coupling between a photon and a diquark with flavour structure ab is in 
impulse approximation described by two very similar diagrams, each with the 
photon coupled to one of the two quarks, depicted in Fig. ^ The diagram where 
the photon couples to quark a of the scalar diquark is given by the loop integral 



(2^ 



Tr 



Sa{k++)r^{k++, k+-)Saik^-) 



xrD,ab{k+--,-k_)CSb{k^)CrD,ba{-k-,k^+) 



(6) 



where Q is the photon momentum, and the incoming and outgoing diquarks 
have momentum P ^ Q/2; the internal momenta are k- = k — P/2, and = 
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Figure 1. Triangle diagrams representing the diquark-photon interaction in impulse approxi- 
mation. 



k + P/2 — Q/2\ the (repeated) flavour index after the semicolon in I^'''"'{P, Q) 
indicates which of the quarks is coupled to the photon. The expression for the 
other diagram is very similar. These two contributions have to be added together 
with the appropriate quark charge. Notice the difference with the expression for 
the meson-photon coupling: there is no factor Nc multiplying the loop integral. 
However, this factor Nc is also absent from the normalisation condition for the 
diquarks (when compared to that of meson bound states), and thus effectively 
this difference cancels out; the calculation for the scalar diquarks and that for 
the pseudoscalar mesons is essentially the same. 

We can define a form factor associated with each of these two diagrams 

2P,FaUQ^) = lf'''{P,Q). (7) 

Thus the electromagnetic form factor for a diquark with flavour structure ab is 
given by 

FabiQ^) = QaFab;a + QbFat;b, (8) 

and the corresponding charge radius is deflned as 



F{Q^) 



(9) 



For the quark propagators and diquark BSAs we use the solutions of the 
rainbow DSE and ladder BSE discussed above. Current conservation demands 
that also the quark-photon vertex F^ is dressed. In impulse approximation, cur- 
rent conservation is guaranteed as long as the quark-photon vertex satisfies the 
vector Ward-Takahashi identity jJTj. The solution of the inhomogeneous ladder 
BSE 

^ ^ G{k^) D'^^%k)^p S{q+) rp{q+, -q.) S{q.) 7. , (10) 
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with the same effective interaction as used in the rainbow DSE for the quark 
propagator does satisfy this constraint, independent of the details of the effective 
interaction TO. 

3 Numerical results 

For the effective quark- antiquark interaction, we employ the Ansatz ^HI 



fc2 



iln 



with J^(s) = (1 — exp 7m = 12/(33 — 2A^j), and fixed parameters 

nit = 0.5 GeV, r = e^ - 1, Nf = i, and Aqcd = 0.234 GeV. The remaining pa- 
rameters Lo and D were fitted in Ref. 15 to reproduce /^r and the chiral 
condensate; however, it turns out that the light mesons are only sensitive to 
the combination Div. The current quark masses rriu = nifi = 3.7 MeV and 
rus = 85 MeV at /x = 19 GeV were fitted to the pion and kaon mass respectively. 
Using the one-loop expression to scale these masses down to = 1 GeV gives 
mu,d{l GeV) = 5.5 MeV and m,(l GeV) = 125 MeV. 

The ultraviolet behaviour of this effective interaction is chosen to be that of 
the QCD running coupling a(A;2); the ladder-rainbow truncation then generates 
the correct perturbative QCD structure of the set of DSEs. In the infrared region, 
the interaction is sufficiently strong to produce a realistic value for the chiral 
condensate of about (240 GeV)^. 

3.1 Diquark masses 

With this model, we have solved the BSE for scalar diquarks. In addition to 
the diquark masses, we also calculated the axial-vector projection of the diquark 
wave function at the origin, which is the analogue of the leptonic decay constant 
of the mesons. To be specific, we define an axial-vector colour and fiavour anti- 
triplet operator, and calculate 



fD,ab = Z2 Tr 



A 
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(12) 



This quantity fD,ab is a purely theoretical object, and does not correspond to 
any physical decay constant of the diquarks. However, it is renormalisation- 
point-independent, just like the meson decay constants. The difference with the 
expression for the meson decay constants is the absence of a factor of Nc- Hence, 
with this definition, we expect fD,ab to be about \/3 times smaller than the 
corresponding pion and kaon decay constants (remember that this factor of Nc is 
also absent from the diquark normalisation condition). Any deviations from this 
expectation would signal dynamical differences between the relative momentum 
behaviour of the diquark and meson BSAs. Our results are given in Tabled 

The light meson properties seem to be rather independent of the details 
of the effective interaction, as long as the interaction generates the observed 
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Table 1. Overview of our results for the pseudoscalar meson and scalar diquark masses and 
decay constants for different model parameters. All quantities except D are in GeV. Numerical 
errors are of the order of 1% for the mesons, and 2% to 5% for the diquarks (unless noted 
explicitly), with the heavier states having the larger numerical uncertainty. 



input parameters 




meson 






diquark 




Lu D [GcV^] 






rriK 


fK 


rriud 




mqs 




0.30 1.25 


0.138 


0.1305 


0.494 


0.1537 


0.988 


0.124 






0.40 0.93 


0.138 


0.1305 


0.495 


0.1537 


0.821 


0.128 


1.10 


0.17(2) 


0.42 0.888 


0.138 


0.1308 


0.496 


0.1541 


0.791 


0.129 


1.07 


0.170 


0.45 0.829 


0.139 


0.1306 


0.496 


0.1544 


0.744 


0.128 


1.02 


0.163 


0.50 0.79 


0.138 


0.1307 


0.495 


0.1558 


0.688 


0.127 


0.96 


0.158 



amount of chiral symmetry breaking JHI- This can be achieved by keeping 
Dijj ~ (0.72 GeV)'^ fixed: the hght meson masses (and other calculated prop- 
erties) remain unchanged for 0.3 GeV < u; < 0.5 GeV, provided D is adjusted 
accordingly, see Tabled However, the diquark masses are more sensitive to de- 
tails of the effective interaction: a 10% increase in lv while keeping Duj constant 
reduces the scalar diquark mass by about 10%. Furthermore, if w becomes smaller 
than about 0.38 GeV (and D larger than about 1 GeV^), we can no longer find 
a stable us or ds diquark. 

Also some of the heavier mesons are sensitive to the range of the strong 
infrared attraction in the model. However, this dependence on u) is not uniform. 
For example, changing to from 0.4 GeV to 0.5 GeV reduces the mass of the scalar 
meson by 10% In contrast, the mass of the first radial excitation of the pion is 
increased by 14% under a similar change The vector mesons are insensitive 
to these changes in lo and D jl5j . 

A complication with these calculations is that the rainbow truncation of the 
quark DSE in general leads to a pair of complex conjugate singularities in the 
quark propagator |191 120j . These singularities are most likely artefacts of the 
model, and progress is being made in eliminating these artefacts |201l21j . Within 
the model under consideration here, the exact location of these singularities in 
the complex plane depends on the parameters ;2ri , but typically the integration 
domain spanned by the BSE encroaches on these singularities for bound state 
masses of slightly above 1 GeV. This limits the exploration of heavier states, and 
it also limits the range of for which we can reliably calculate the electromag- 
netic form factor. 

3.2 Electromagnetic form factors 

Next, we use the obtained diquark BSAs to calculate the electromagnetic form 
factors. For completeness, we also include the charge radii for the pseudoscalar 
mesons, illustrating that the pion and kaon properties are indeed insensitive to 
the details of the interaction. For the various mesons we have 
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FkAQ^) = \fu-sAQ'') + \fu-s--s{Q'') ^ (14) 

F^o(Q2) = AFd-s;d{Q'') + \Fd-s.-s{Q^), (15) 
whereas for the diquarks we have 

FudiQ'') = \FudAQ^)-^F^dAQ^), (16) 

FusiQ^) = \FusAQ^) -^Fus-AQ^) ^ (17) 

FdsiQ^) = -\FdsAQ^)-\Fds-AQ^)- (18) 

Impulse approximation, in combination with the rainbow-ladder truncation, en- 
sures current conservation: each of the form factors Fqq-g{Q^) and Fgq-g{Q^) is 
one at = 0, independent of the flavour label q. Thus we obtain the correct (in- 
teger) charges for the physical mesons, whereas for the diquarks we get fractional 
charges 1/3 for the ud and us diquark, and —2/3 for the ds diquark. Note how- 
ever that these diquark states are coloured objects, and therefore confined, and 
these fractional charges cannot be observed directly. With the above definition 
of the "physical" form factors, i^afe(O) gives the charge of the state under con- 
sideration. This (fractional) charge is divided out in the definition of the charge 
radius, see Eq. Q. 

Table 2. Overview of our results for the pseudoscalar meson and scalar diquark charge radii 
squared, all in fm^ (the flavour label q stands for both u and d quarks). The estimate of the 
combined numerical errors is less than 0.01 fm'^ for the meson charge radii, about 0.01 fm^ for 
the ud radii, and about 5% to 10% for the us and ds radii. For lo = 0.4 GeV, D = 0.93 GeV^, nu- 
merical difhculties due to the complex singularities in the quark propagator prevent an accurate 
determination of the charge radius for the us and ds diquarks. 



input 








meson 










diquark 








D 


' TT 


qs\q 


r^- - 

' qs;s 








^2 

qa;q 


' qs\s 


' US 




0.40 


0.93 


0.433 


0.460 


0.209 


0.376 


-.084 


0.51 










0.42 


0.888 


0.433 


0.458 


0.210 


0.375 


-.083 


0.51 


0.45 


0.27 


0.63 


0.36 


0.45 


0.829 


0.438 


0.460 


0.211 


0.377 


-.083 


0.51 


0.49 


0.27 


0.71 


0.38 


0.50 


0.79 


0.441 


0.455 


0.211 


0.374 


-.081 


0.51 


0.49 


0.27 


0.71 


0.38 



Our numerical results for the charge radii are given in Table [2 We work in 
the isospin symmetry limit, and thus we do not discriminate between the u and 
the d quark, except for the obvious difference in their electromagnetic charge. As 
already mentioned, the pion and kaon charge radii are insensitive to variations 
of the parameters as long as the integrated infrared strength is kept fixed. Also 
the diquark charge radii turn out to be relatively insensitive to these variations, 
despite the fact that the diquark masses change considerably over this range in 
the parameters. This is at least partly due to the fact that the vector meson 
masses are also insensitive to these variations in the parameters, and they show 
up as poles in the dressed quark-photon vertex. 



8 



Electromagnetic properties of diquarks 




photon [GeV^] photon q' [Gcv'] 



Figure 2. The pion (left) and ud diquark (right) form factor, together with the fitted forms, 
for CO = 0.42 GeV, D = 0.888 GeV^ The analytic forms are fitted in the range -0.3 GeV^ < 
< 1.5 GeV^. 



Indeed, all of the form factors do show a pole-like structure in the timelike 
region, see Fig. |2 For the pion, a monopole fit to the form factor 

using 27 calculated values of F{Q'^) with in the range -0.4 GeV^ < Q2 < 
1.6 GeV^ gives a mass-scale that is very close to the p meson mass in this model, 
m?p = 0.55 GeV^, with a of about 0.001. Again, these results are almost 
independent of the details of the interaction, as can be seen from Table 01 

For mesons (and diquarks) involving one strange quark and one up or down 
quark, one has to choose whether to attempt a monopole fit to the form fac- 
tors corresponding to the individual diagrams of Fig. ^ or to the physical form 
factors. The neutral kaon form factor cannot be fitted by a simple monopole; 
however, it might be fit by a difference of two monopoles. Indeed, the vector 
meson dominance [VMD] mechanism suggests that one should fit the individ- 
ual diagrams of Fig. ^ it is the vector meson pole in the quark-photon vertex 
that is at the origin of the success of VMD models. Thus, for the kaon, we have 
fitted Fus;u{Q'^) and Fus;s{Q'^) over the same range, and the agreement with a 
monopole is reasonable, though not as good as for the pion; the value is about 
0.01. Furthermore, the masses deviate somewhat more from the vector mesons 
TTip = 0.55 GeV^ and = 1.15 GeV^, see Table 01 As the pole positions are 
approached, these form factors become increasingly better represented as VMD 
monopoles with the p and 4> masses jl61 122| . 

For the diquarks, a monopole fit is not very good, as can be seen from the 
right panel of Fig. [2 the of the monopole fit is around 0.1 for the ud diquark. 
A much better fit can be achieved by using the two-parameter fit form 

which gives a of 0.001. The preferred power appears to be a = 1.4, independent 
of the details of the interaction. However, this does depend somewhat on the 
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Table 3. Estimates for monopole fit masses in GeV^ (the flavour label q stands for both 
u and d quarks). 



input parameters 




meson 




diquark 




ud diquark best fit 


uj [GeV] 


D [GeV^] 


pion 


qs;q 


qs;q 


ud 


qs;q 


qs; s 


mass 


power 


0.40 


0.93 


0.545 


0.529 


l.fl6 


0.51 






0.62 


f.35 


0.42 


0.888 


0.545 


0.530 


l.f09 


0.5f 


0.4 


0.8 


0.62 


f.35 


0.45 


0.829 


0.540 


0.527 


l.fOO 


0.50 


0.54 


0.85 


0.63 


f.38 


0.50 


0.79 


0.538 


0.530 


l.f03 


0.50 


0.5f 


0.89 


0.62 


f.37 



range that is used in the fitting. Close to the p pole in the timelike region, it does 
approach a pure VMD monopole form. A similar fit to the calculated pion form 
factor gives a power of 1.02, which is hardly an improvement over the simple 
monopole fit. These fits with an arbitrary power can only be good for low and 
intermediate values of Q^, since we know that asymptotically, the pseudoscalar 
meson and scalar diquark form factors vanish like 

For the qs;q and qs;s form factors, with q standing for u or d quarks, the 
monopole masses in TablelHlare an estimate based on the low-Q^ behaviour only, 
rather than a fit over the same range of Q^. Complications with the analytic 
continuation of the quark propagator solution necessary for these calculations, 
and in particular the complex singularities in the quark propagator, prevent us 
from obtaining the form factor for these heavier states away from ~ 0. 

In conclusion, we can say that simple VMD models work very well for the 
pion form factor, but they are less and less accurate for the heavier mesons. Nev- 
ertheless, the monopole forms for the diquarks are a reasonable approximation to 
these microscopic form factor calculations, and might be useful for applications 
in baryon calculations. An interesting consequence of these monopole fits is that 
the us diquark form factor could potentially have a zero crossing in the spacelike 
region: if the mass-scale in the Fus;s fit is more than twice as large as the mass 
scale in the Fus-u fit, the "physical" form factor ^Fus-^uiQ"^) — ^Fus-siQ'^) will 
approach zero from below for ^ oo. (In terms of VMD masses that means 
if M| > 2M^.) However, both VMD and our current calculations indicate that 
this is not the case for realistic up, down, and strange quark masses. 

4 Concluding remarks 

The calculated charge radius of the ud scalar diquark, = 0.71 fm, is 8% larger 
than the pion charge radius, r^^ = 0.66 fm. This suggests that these diquarks are 
about 8% larger in size than pions, which agrees with the observation that the 
diquark BSA is about the same amount narrower in momentum space than the 
pion BSA These results appear to be insensitive to the details of the effective 
quark-quark interaction kernel, even though the actual diquark masses are sensi- 
tive to details of the interaction. Our results disagree with calculations within the 
Nambu-Jona-Lasinio [N JL] model [23 , which suggested that the electromagnetic 
form factors of the pion and ud scalar diquark are almost identical. Furthermore, 
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the obtained charge radius in those NJL calculations is significantly smaller than 
the measured pion charge radius. 

Our result for the ud scalar diquark charge radius is about 20% smaller than 
the experimental proton charge radius, Tproton = 0.87 fm. If we take these radii as 
an indication of their physical size, it is indeed quite plausible to picture baryons 
as bound states of a quark and a diquark. For quantitatively reliable calculations 
of electromagnetic properties of nucleons however, we need not only the scalar 
diquark, but also the axial-vector diquark properties. 

Recently, it has been suggested that pentaquarks might be composed of two 
diquarks (coupled together in a colour triplet) coupled to an antiquark An- 
other possibility involving diquarks would be a (colour anti-triplet) diquark cou- 
pled to a triquark {qqq) in a colour triplet state [2]. One would expect pen- 
taquarks to be somewhat larger in size than the ordinary baryons. Based on 
the charge radius of the diquarks and that of the proton, we would estimate 
the radius of a pentaquark (in a diquark-diquark-antiquark configuration) to be 
''pentaquark ~ 1 fm, or about 15% larger than that of the proton. 

In the timelike region, all form factors can be approximated quite well by 
a VMD like monopole. For increasing spacelike values oi , a VMD monopole 
becomes less and less accurate, in particular for heavier mesons and diquarks. For 
moderate spacelike values of Q^, the ud diquark form factor may be interpolated 
by [M2/(g2 _^ m2)]1-4, with a mass of about M ^ 0.79 GeV. Asymptotically 
however, we expect this form factor to vanish like just like the pion form 

factor. 
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